Order reduction in semiclassical cosmology. 



Grazyna Siemieniec-Oziehlo and Andrzej Woszczyna 
Astronomical Observatory, Jagellonian University, 
ul. Oria 171, 30-244 Krakow, Poland 
woszcz@oa. uj. edu.pl 

We investigate the Robertson- Walker cosmology with Lagrangian 
where L ra d means classical source with traceless energy-momentum tensor. We weaken the self- 
consistence condition (L. Parker, J. Z. Simon, Phys. Rev. D47(1993),1339). Quantum corrections 
are expressed as contributions to the effective equation of state. We show that the empty space-time 
is stable within the class of radiation-filled expanding universes with no order reduction of the field 
equations. 

PACS numbers: 9880 



I. INTRODUCTION 

Controversial fourth order differential equations, which govern the semiclassical cosmology can be reduced to second- 
order [jlj , [^| , and in this way, exempted from quantum-originated instabilities || [Q . The reduction is based on the 
self-consistence condition, i.e. the assumption that both equations and solutions arc pcrturbatively expandable in 
h. Under this condition the universe becomes an ordinary mechanical system with a two-dimensional phase-space 
corresponding to the single degree of mechanical freedom - the scale factor a(t,Ti). Self-consistent theory is still 
renormalizable H, Minkowski space-time regains stability in the class of homogeneous and isotropic models, quasi- 
inflationary phenomena disappear Similar reduction techniques are being applied to gravity with higher than 
fourth-order derivatives |7| and also in other branches of physics [0 . 

However, imposing the self-consistence condition on the cosmological scale a(t) encounters some difficulties. In 
a universe with vanishing spatial curvature still remains a freedom to multiply metrics by an arbitrary constant 
factor, therefore the scale a(t) is not a measurable quantity. Requirement for a(t, H) to be ^-expandable is physically 
unclear. In open or closed universes this freedom is reduced by the choice k = ±1, and the scale factors are uniquely 
determined by cosmological observables: the Hubble parameter H and the energy density e. Yet arbitrarily small 
changes in any of these observables in the vicinity of critical density e — A = 3H 2 may result in indefinite changes of 
a(t). The perturbative character of the energy-momentum tensor (and consequently the equations) with respect to 
an arbitrary chosen parameter, in general, does not imply the same property of the metrics. Finally, expanding a(t) 
in the equations, which contain fixed curvature index k limits quantum corrections to only those, which preserve 
the same sign of the space curvature. This limitation is particularly severe for a flat universe, where generic quantum 
corrections would contribute to the space curvature unless the k — condition prevents that. This limitation cannot 
be derived directly from the Lagrangian and, in fact, it forms an additional constraint imposed on the theory (which 
is not even true in classical gravity^). 

Not arguing with the very idea of self-consistence, we draw attention to some circumstances which are important 
for semiclassical cosmology: 

1. Without harm to the reduction procedure, one can release the consistence condition for the scale factor, de- 
manding instead the same property for cosmological observables (the Hubble parameter, etc.) 

2. For a radiation filled universe with vanishing cosmological constant A = the self-consistence condition is 
superfluous, since the original equation is of second (!) order. Terms with higher derivatives cited by classical papers 
Q contain an additional hidden factor h and are eventually eliminated in the first order expansion. 

3. We show that quantum corrections form the equation of state of a barotropic fluid, and discuss the stability of 
the Minkowski space-time on the ground of dynamic systems theory. 



1 Note that the Lemaitre universes, which are of positive space curvature, are obtained from flat universes (not closed!) when 
A diverges from zero. 
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II. CONDITION OF SELF-CONSISTENCE FOR HUBBLE'S EXPANSION RATE. 



We consider semiclassical gravity theory with the Lagrangian R+aihR 2 +a2hR^R^ l/ + L ra d, where L ra d represents 
classical radiation or another thermalizcd field of massless particles. Typically, cosmologies containing the R 2 and the 
WRuv terms lead to 4 order equations || and violate the stability of empty space ||. 

We write quantum terms on the right-hand side of the field equations and treat them as corrections to the energy- 
momentum tensor. We think of fi as the theory parameter, which can take arbitrary values, so the limit transition 
7i — ^ defines the classical limit of the theory. The field equations we write in the Einsteinian form R^ v — \Rg^ v + 
ApjU!/ = but with the modified, effective energy-momentum tensor 
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and the constant a 3 is some combination of a\ and 0:2 (Robertson- Walker symmetry have been partially exploited to 
derive formula (Q). For more precise explanation see ||.) Derived in this way the (0,0)-equation 
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contains four fundamental constants, two of them classical - the gravitation constant k (further on we put 8ttk = 1), 
the cosmological constant A, and two quantum ones - a\ and 0:3. There are also two quantities which define a 
particular solution: the constant of motion /j, = £o a oi an d the index of space curvature k. Therefore, the transition 
from classical to quantum theory with the self-consistence a — ao + Tia\ imposed on (^) preserves the type of space 
curvature, including the strong k = limitation for the flat universe. 

One can get rid of the last two constants, and consequently, of the constraints they bring, by introducing the Hubble 
expansion parameter H = . Differentiating (Q) twice, we obtain the fourth order equation for H, which contains 
only fundamental constants^] 
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(5) 



This equation describes the dynamics of Robertson- Walker models with arbitrary space curvature, and what is equally 
important, it is expressed in terms of observable quantities. A self-consistence condition imposed on measurable 
quantities has well defined physical meaning. We adopt Simon's ansatz to H, namely we state that H(t) — H c i ass (t) + 



2 The reverse procedure would give the equation with two parameters of continuous values. Consequently, the equation (g) 
formally has a broader class of solutions than (^). However, the freedom to choose k as different from 0, ±1 is a trivial one, 
and resolves itself to rescaling the metrics by a constant factor. 
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hH quant (t) is perturbative in h. Now, the procedure of the order reduction can be done in two ways: 

1) one can differentiate twice the zeroth-order expansion (equation (|^) with at = a 3 = 0) to find the third and fourth 
derivatives and eliminate them from the full equation (|5|) - this is equivalent to what is done in 0], 

2) substitute the expansion H(t) = H c i ass (t) + hH quant (t) directly into (|E|) and abandon terms second order in h or 
higher . 

In both cases we obtain the second order equation 

d 2 H dH , , 
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dt v ' 



dt 2 
d 2 H 



dt 2 



dH 



3 (51fiai + fia 3 ) ^f- + 3 (90?iai + Tia 3 ) H 2 - 2A (67kxi + ha 3 ) 
dt 



AH 



A59h ai 



dH 



dH 



v 3 (372^ai - ha 3 ) H 2 - 2A (697iai + Has)) 
A (3 (180toi - ha 3 ) H 4 - 20AKh ai H 2 + 8A 2 hai) , 



(6) 



which is nonlinear both in H and its derivatives. So strong nonlinearity allows one to find exact solutions only in some 
particular situations. The is not the case in equation (^). However, this equation becomes much more transparent 
after one rewrites the quantum corrections as contributions to energy density and pressure. Qualitative analysis is 
then enabled. 

Let e and P denote respectively effective energy density and effective pressure, i.e. each of these quantities is 
supplemented by quantum corrections. The universe dynamics is determined by the system of the Raychaudhuri (Q) 
and the continuity (@) equations 



dH 
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(h_ 

dt 



= -H 2 - - (3P + e - 2A) 
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(7) 
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We differentiate (0), substitute into (g) and apply the continuity equation (g) to get the relation between pressure, 
energy and the cosmological constant in differential form 
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As a matter of fact, one can solve equation (||) analytically, however the solution takes unclear implicit form. This is 
much simpler to follow the other way. The solution of (|^) must be a function of the energy density and cosmological 
constant solely, hence P(e, A) is independent of the expansion rate H. Therefore the limit transition H 2 — > 3(e + A) 
does not affect its values, and the general solution is identical with the integral found for the flat universe. In the last 
case the equation 
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is an analogue to equation (^). Its order reduces by two, and finally the equation takes a particularly simple form 

dH 2 2a 3 h 



dt 
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Now, comparing (^) with the Raychaudhuri equation (Q) we obtain the equation of state of cosmological substratum 
in the form of the algebraic relation 



P=-e 
3 



Aa 3 h 



A 



(12) 



Function Pfe, A), defined by (12), fulfills the differential equation(|9j) with an accuracy to terms o(h). By simple 
calculation , one can confirm that the exact solutions found by Parker and Simon also obey ([lj) . 

As we have already mentioned, the equation of state ([l2]) is barotropic, i.e. effective pressure is solely the function 
of the effective energy density (including the energy of vacuum A). While reducing the equations order we eliminate 
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contributions to the energy- momentum tensor coming from the expansion rate ||; therefore the universe evolution 
becomes a reversible process (equations (@)-(|8|) are invariant under the time reflection t — > —t). 

Quantum corrections contained in Jl2[), and consequently the dynamical system are ^ ree °f the a i constant. 
The only term multiplied by a\ which survives the reduction procedure H) , has been assimilated here by the effective 
energy density!]. 



III. THE A = CASE. 



Its worth noticing that in some physically interesting situations the reduction procedure eliminating higher order 
derivatives is redundant. In the radiation filled universe with null cosmological constant the correction haiH^]) , which 
formally appears as linear in h, actually is quadratic, and consequently should be abandoned as the o{K) term. To 
show this let us express the traceless tensor ^Hjf in terms of the Ricci scalar and the effective energy density 

The field equations with the energy-momentum tensor (|l|) show that the scalar R involves the trace of the tensor 
WHp, namely R = a 5 h (5) H£, so it is a quantity linearly dependent on h. 

Writing ( 3 'H£ in terms of the effective energy density e with the accuracy to terms o(T%) we get R = ^a^he 2 . Tensors 
^'Hf} and ^Hj} can be rewritten as 
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Now it is clear that only the second of the expressions aifr^Hp and a^fS 3 ^ H£ is essentially linear in % and forms 
the first-order quantum contribution to the energy-momentum tensor. The first one aifir'Hf}, which carries all 
higher derivatives is actually square in h. This is closely related to the absence of particle creation in the radiation- 
filled Robertson- Walker universe (see ]TT[] and papers cited there.) The theory with the energy-momentum tensor 



m rp(rad) 



leads to the effective equation of state P = he 



4a 3 h r 2 



e , which is perfectly consistent with 



IV. STABILITY OF THE EMPTY SPACE - DYNAMICAL SYSTEMS APPROACH. 



The equation of state of the form P = P(e, A) (or more generally P = P(e, A, H) sec |l2[ ) uniquely determines cos- 
mological evolution. The system (0^), which defines the universe dynamics in the {H, e}-phase space is autonomous. 
Choosing a point in the {H, e}-phase space, one determines uniquely the metrics in the initial moment as well as the 
metrics' evolution in time.^j The stability of the Minkowski space-time is defined by the stability of the {H, e) = (0, 0) 
point in the {H, e}-phase space under the condition A = 0. For the equation of state ( |l2|) discussed in the preceding 
section the autonomous system (0-||) reads: 



3 In this approach the quantum corrections modify effective energy density and pressure, not the fundamental constants like 
in | 

4 We abandon here a trivial freedom to multiply the flat universe metrics by a factor constant in time. 



4 



dH 
~dt 
de 
~dt 



= -H 2 - 



— = -4H 



1 % 2 

3 J 



and its trajectories form levels of the integral 
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(13) 
(14) 

(15) 



The phase portrait of the system (|13[|1J) is shown on Fig. 1. For completeness and also for readers convenience, we 
attach Fig. 2. showing classical Friedmanian dynamics in the same representation. The phase structure of classical 
radiation-filled universes and the phase structure defined by ( fu^li!) are topologically equivalent in the low energy 
limit. This is so because one cannot enrich the structure of the {H, e}-phase plane without violating the standard 
energy conditions. On the other hand, according to ( p"2| ) these conditions are well fulfilled for low and positive energy 
densities. The equation of state (|l^) formally admits violation of the energy conditions but these states appears 
already in the Planckian regime and hence, far beyond the region where semiclassical approximation is valid. (The 
dotted region in the upper part of Fig. 1, which contains three 'additional' critical points must be recognized as 
nonphysical) . 









if 












FIG. 1. Dynamics of semiclassical radiation 
filled universes. 



FIG. 2. Dynamics of Friedmanian (classi- 
cal) radiation filled universes. 



An essential property of the system (|l3-14) is the absence of solutions that change the energy density from positive 
to negative, or the reverse. (Such behaviour was possible in the original semiclassical theory and disqualified the 
empty space as a ground state.) Indeed, on the strength of ([jjj), the initial condition e = results in e + P = 0, 
and consequently the right-hand side of equation ( |l4| ) vanishes. Both conditions e = and de/dt = ensure that the 
state of the zero energy density is 'persistent'. This is consistent with the results based on the functional integral 
formalism where all higher derivative terms responsible for instability are eliminated by regularisation of the 
energy-momentum tensor. 

The stability of Minkowski space-time is the same as in the classical theory. In both cases, the classical or the 
quantum, the (H, e) — (0, 0) point is a three-fold point with elliptical sector jyj and its type does not depend on 
the value of h. This means that the phase space is structuraly stable against quantum corrections in the low energy 
density limit. This nontrivial property does not follow from the solutions analicity in h, but from the form of the 
energy density tensor (|l|) f\ 



3 In general, a three-fold point may bifurcate into simple critical points under smooth changes of the equation coefficients. 
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V. SUMMARY AND CONCLUSIONS 



In the reduced Simon-Parker theory the energy-momentum tensor is renormalized to take the hydrodynamic form 
with a simple, barotropic equation of state. 

The self-consistence conditions for semiclassical cosmology can be imposed on observable quantities and weakened. 
By demanding the Hubble expansion rate to be perturbative in Ti we allow the space curvature to alter from while 
quantum corrections to the flat universes occur. 

In the particular case of the radiation-filled universe and vanishing cosmological constant, the dynamics of the 
Robertson- Walker universe in the (original) semiclassical theory is described by a second order equation, therefore 
it does not need either the reduction or additional conditions of self-consistence. The reason lies in the absence of 
particle creation in the radiation filled universe, which manifests itself as an additional factor fi 'hidden ' in the tensor 
^Hfn,. This eventually eliminates all the higher derivative terms. 

Minkowski space-time has the same stability character as for Einsteinian gravity, which is consistent with results 
based on the functional integral formalism Jl3| ] . The stability of Minkowski space-time is independent of the numerical 
value of the Planck constant. In the language of dynamical systems theory, this property is called the structural 
stability of the {H, e}-phase space against changes of Ti. 

Its worth noticing that the Liapunov stability of the environment with equation of state (|l2| ) with respect to 
position-dependent perturbations is also the same as for the classical radiation-filled universe ]16[, in contrast to the 
original semiclassical theory, where quantum corrections let inhomogeneities grow. This suggests an insignificant role 
for semiclassical corrections in the processes of structure formation in the early universe. 
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